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Einstein-Podolsky-Rosen (EPR) steering is a type of quantum correlation which allows one to remotely pre¬ 
pare, or steer, the state of a distant quantum system. While EPR steering can be thought of as a purely spatial 
correlation there does exist a temporal analogue, in the form of single-system temporal steering. However, a 
precise quantification of such temporal steering has been lacking. Here we show that it can be measured, via 
semidefinite programming, with a temporal steerable weight , in direct analogy to the recently proposed EPR 
steerable weight. We find a useful property of the temporal steerable weight in that it is a non-increasing function 
under completely-positive trace-preserving maps and can be used to define a sufficient and practical measure of 
strong non-Markovianity. 
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Quantum entanglement, Einstein-Podolsky-Rosen (EPR) 
steering, and Bell non-locality are three of the most intrigu¬ 
ing phenomena in quantum physics and, in varying degrees, 
are thought to act as resources; fuel that powers a range of 
quantum technologies. Entanglement [1-3] comes in hand- 
in-hand with the complexity of quantum systems, and may be 
behind the potential speed-up of quantum computation. Bell 
non-locality and EPR steering are thought to be the driving 
power of quantum cryptography, and have both been recast 
in that language. For example, in a quantum key distribu¬ 
tion scenario, two parties wish to generate a secret key using 
shared quantum states as a resource. If one party (Bob) trusts 
his own experimental apparatus but not that of the other party 
(Alice), a violation of a steering inequality [2^f] can be used 
to certify that true quantum correlations exist between their 
shared states. In stricter terms, such a test proves to Bob that 
the correlations he observes between his measurement results 
and Alice’s cannot be described by a local hidden state model; 
his state is truly being influenced by Alice’s measurements in 
a non-local manner. As with entanglement, one quantify the 
amount of steering that is possible with a given shared state 
via a range of possible measures [5-8]. Very recently, a pow¬ 
erful example of such a measure, the steerable weight, was 
proposed by Skrzypczyk et al. [9, 10]. 

In EPR steering the notion of non-locality, via space-like 
separations between parties, plays an important role. If we re¬ 
lax this constraint, and consider time-like separation of mea¬ 
surement events, can the concept of steering still be used in a 
meaningful way? We can find inspiration in the fact that there 
do already exist other types of non-trivial temporal quantum 
correlations complementary to both Bell non-locality and en¬ 
tanglement. For the former, one of the most well-known ex¬ 
amples is the Leggett-Garg (LG) inequality [11], which can 
be used to test the assumption of “macroscopic realism”, in 
contrast to the non-local realism tested by Bell’s Inequality, 
and for which experimental violations have been observed in 
a large range of systems [12-14]. For the latter, motivated by 


the Choi-Jamiolkowski (CJ) isomorphism [15], which equates 
the correlations in a bi-partite quantum system with two-time 
correlations of a single quantum system, the notion of tempo¬ 
ral entanglement has been proposed in various forms [16-22]. 
Returning to steering, the concept of temporal steering, and a 
temporal steering inequality, was recently introduced by Chen 
et al. [23]. Also inspired by the CJ isomorphism, they showed 
that, even without the assumption of non-locality, the concept 
of one party not trusting the earlier measurements made by an¬ 
other party delineates between certain classical and quantum 
correlations. Not only does this have direct practical applica¬ 
tions in verifying a quantum channel for quantum key distri¬ 
bution (QKD), it was recently shown that temporal steering, 
like EPR steering [24, 25], is intimately linked to the concepts 
of realism and joint measurability [26-29]. 

Still lacking however is a measure to quantify these “tem¬ 
poral steering” quantum correlations. Here, in analogy to 
the EPR steerable weight [9, 10], we define the temporal- 
steerable-weight (temporal-SW) as a measure of temporal 
steering. We prove that the temporal-SW is non-increasing 
under a completely-positive trace-preserving (CPT) map and 
can be used to define a sufficient but not necessary measure 
of non-Markovianity. In the same way that the spatial steer¬ 
able weight can be considered a measure of strong entangle¬ 
ment, since not every entangled state is steerable, we define 
the temporal-SW as a measure of strong non-Markovianity 
because it vanishes for weak non-Markovian process. (This 
is also in analogy to, e.g., the phenomenon of strong non- 
classicality, which can be detected and quantified by a weaker 
criterion of non-classicality [30, 31]). We show this by com¬ 
paring the non-Markovianity measured by the temporal-SW 
to an existing entanglement-based measure [32], and find that 
it is, as expected, less-sensitive. However, the temporal-SW 
is, in principle, easier to implement experimentally, as it does 
not require the use of an ancilla, nor full process tomography. 
These results, together with with a few illustrative examples 
discussed mainly in the Supplementary Material [33], suggest 
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FIG. 1. (Color online) Schematic diagram of temporal steering. In 
the beginning, Alice performs the measurement F a \ x = mLm „\ x 
on an initial state po . Then, po is mapped to p a \ x and sent into a 
quantum channel A. Finally, Bob receives the assemblage {cr^ x } at 
time t. 


that temporal steering can serve as a unique and useful quan¬ 
tum resource. 

Temporal steerable weight .— Now we introduce the con¬ 
cept of a temporal steerable weight in analogy to the spa¬ 
tial steerable weight introduced recently by Skrzypczyk et 
al. [9, 10]. In the standard (spatial) EPR steering scenario, 
Alice performs a POVM (positive-operator valued measure) 
measurement F a \ x = M^ x M a \ x , ]T a K\x = 1, on a state 
Pab shared with Bob and creates the assemblage {a a \ x }, 
where a is the measurement result and x is the basis of the 
measurement. In defining a steering inequality, or a steerable 
weight, one assumes that Bob does not trust Alice, nor her ex¬ 
perimental apparatus, and wishes to distinguish between true 
manipulation of his local state via quantum correlations and 
correlations that cannot be distinguished from some classi¬ 
cal theory, typically a local hidden state model. In tempo¬ 
ral steering, we also let Alice perform a POVM measurement 
F a \ x = M^ x M a \ x but on a single system in an initial state 
po at time t = 0. After the measurement, the initial state is 
mapped to p a \ x (see Fig. 1): 


PO '“A Pa\x 


M a \x PO ^ a \ x 
p(a\x) 


( 1 ) 


with the probability p(a\x) = ti(M a \ x poM^ x ). After this 
initial measurement, the state p a \ x is sent into a quantum 
channel A for a time t. At time t. Bob receives the system 
and performs quantum tomography to obtain the state er a | x , 
i.e., A(p a | x ) = cr a | x . To mimic the unnormalized assem¬ 
blage [9, 10] in standard EPR-steering, we define the unnor¬ 
malized states in temporal steering 

a l\x = P( a \ x ) a a\x, (2) 

where the superscript T reminds one that the assemblage 
{<cr^, } is for temporal steering. 

However, the quantum channel may be noisy, obliterat¬ 
ing the influence of Alice’s measurement choice, or Alice’s 
measurement results could have been fabricated via classical 
strategies. In these cases, o^ x may include, or be entirely 
described by, an unsteerable assemblage which we define as 


= E p ( A ) p («i*i A ) (3) 

A 


where P(A) = 1. We have written the result a, condi¬ 
tional on the basis x, with a subscript notation a\ x = a\x. In 
the EPR setting A represents a local hidden variable which de¬ 
termines the possible correlations between Alice’s and Bob’s 
measurement results from a source which obeys classical re¬ 
alism. As in that case, when Alice reveals her measure¬ 
ment results, Bob can update his knowledge of his state, as 
indicated by two equal forms (by applying the chain rule), 
J2\ p W p (a\ x \X)ax = Y^\ P ( a \ x ) p (Ma\x)v\- Then the 
unsteerable states are those states which obey the classical 
(realism) chain rule for Alice’s joint measurement results, as 
shown in a recent work on steering witnesses [26]. No mat¬ 
ter what happens during the transmission, Bob’s task is to 
check whether the assemblage he receives can be written in 
the hidden-state form [Eq. (3)] or not. If he can, this means 
the state Bob receives is independent of the basis x Alice 
chooses to measure in. As mentioned above, this may be be¬ 
cause the quantum channel is too noisy, such that the influence 
of Alice’s measurements is no longer discernable, or Alice’s 
measurement results could have been fabricated via classical 
strategies. On the other hand, if the assemblage Bob receives 
cannot be written in the form of Eq. (3), he is convinced that 
Alice has influenced his state by her choice of measurement. 
In this case, we call the assemblage Bob receives “temporally 
steerable” and is symbolized as }. 

To determine the steerable weight, one considers the over¬ 
lap between the state Bob receives and the unsteerable assem¬ 
blage, such that his state can be written as a mixture 

(J a\x=P°a\x +( 1 ~^K| :C - W 

To quantify the “steerability in time” for a given assemblage 
{VaLx}’ one has to maximize p, i.e., maximize the proportion 
of crj’|^ s . Then, the “temporal steerable weight” can be defined 
as TSW = 1 — p*, in which p* is the maximum of p and can 
be obtained from semidefinite programming [9, 10, 36]: 

p* = max tr ^ ox 
x 

subject to CqIx — E Dx(a\x)ox > 0 Va,x (5) 

A 

a a > 0 V A, 

where a\ = pax , and D\(a\x) are the extremal determin¬ 
istic values [9] of the conditional probability distributions 
P(a| x |A). Equation (5), which is formulated as a semi-definite 
program (SDP), can be numerically implemented in various 
convex optimization packages, e.g.. Refs. [37, 38]. 

So far, the formalism is parallel to the standard EPR steer¬ 
able weight [9]. The primary difference is that {cr a | x } in 
Ref. [9] is created through the entanglement between Alice 
and Bob. Here, {cr^, } is created through Alice’s measure¬ 
ment and the influence of the quantum channel A. In the Sup¬ 
plementary Material [33], we give an explicit pedagogical ex¬ 
ample of how to evaluate the temporal steerable weight. 

Measure of non-Markovianity. — Now we apply the in¬ 
troduced temporal steering weight as a measure of non- 
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Markovianity. Non-Markovianity is a term used to describe 
the situation when an environment surrounding a quantum 
system has memory of its past evolution. It is an impor¬ 
tant concept both because many natural and man-made quan¬ 
tum systems exist in a regime where the assumption of a 
Markovian (memory-less) environment fails, but also because 
it can lead to counter-intuitive results regarding the decay 
of quantum effects, particularly when the quantum system 
is strongly coupled to the surrounding environment. There 
has been a range of efforts at constructing measures of non- 
Markovianity, typically based on a scenario where the time 
evolution of a quantum system is analysed for non-Markovian 
properties. Arguably, the most popular measures of non- 
Markovianity were introduced in Refs. [32, 39]. Recently, an 
attempt to classify these non-Markovianity measures in a uni¬ 
fied framework was described in Ref. [40]. Useful for us here 
is the approach taken in [39], which is based on observing the 
behavior of the trace distance between two quantum states. 
They derived a measure of non-Markovianity by noting that 
all CPT maps $ are contractions of the trace distance metric, 
and a given dynamic processs is defined as Markovian if the 
map is divisible, i.e. $(r + f, 0) = $(r + 0), for 

all positive t and r. These two properties lead to the mono¬ 
tonicity of the trace distance, and violations of this monotonic¬ 
ity indicate the occurrence of non-Markovian dynamics. In a 
similar way, below we prove that the temporal-SW of a sys¬ 
tem undergoing a CPT map is also a non-increasing function, 
i.e.. 


TSWp > TSW $(t)p (6) 


for a CPT map <!>(t). Together with the property of divisibil¬ 
ity, one can conclude that the temporal-SW decreases mono- 
tonically under Markovian dynamics. Therefore our mea¬ 
sure of non-Markovianity is defined by integrating the positive 
slope of the temporal-SW 


■A/tsw = / dta TS w(t,po,$), 

J <ttsw >0 


(7) 


where (7tsw(C Po, $) = ^TSW$( t ) Po is the rate of change of 
the temporal steerable weight. In the examples discussed in 
the Supplementary Material [33], we demonstrate explicitly 
how one can use this as a practical measure of strong non- 
Markovianity. Here we discuss only the following example. 

Proof of the monotonicity of temporal-SW under Marko¬ 
vian dynamics .— First, we prove that the temporal-SW of a 
system undergoing a CPT map is a non-increasing function, 
as given by Eq. (6). To obtain the temporal-SW of a qubit at 
timefi, one needs the quantity, — D\ 1 (a\x)d\ 1 , 

in which the set {<r Al } is chosen to maximize Tr(J] A cr Al ) at 
time f i. Summing all the measurement outcomes a and taking 


the trace, we have 


Tr 


J2 ^*(*1) (“I®)**! 


a Ai 


( 8 ) 


= Tr 




i * 

- 1 A*1 5 


where we have used the properties ]P a Z) A (a|:r) = 1, and 
Tr ^2 a <£. (ii) = 1. Similarly, to obtain the temporal-SW 
of the qubit at a later time f 2 = ft + r, one also has 


Tr 


D\ 2 (a\x)a\ 2 

_ CL CL A2 


= 1-/4 (9) 


where {oa 2 } is chosen to maximize Tr(^4 d\ 2 ) at time £ 2 - 
One can also perform a CPT map $(r) to Eq. (8), giving 


Tr 


_ a cl Ai 


= Tr 


EE (a|z)($(T)a- Al ) 

. ci CL Ai 


( 10 ) 


Since <1 >(t) is a trace-preserving map, the value of Eq. (10) 
is still 1 — p*. However, we know that the set {<ta 2 } is the 
optimal way to maximize Tr(£4 d\ 2 ) at time t .2 for Eq. (9). 
Therefore, comparing Eq. (9) with Eq. (10) would give 


1 - p\ > 1 - p* 2 , 


(ii) 


This proves the theorem given in Eq. (6). Employing the di¬ 
visibility of Markovian dynamics leads to the monotonicity of 
the temporal-SW: 

TSW{4>(t + t. 0)cr a | x }= TSW{$(r + t, 0K, X } 

<TSW{<b(i,0K|J (12) 

An example of non-Markovianity of a spin-boson 
problem .— Exact solutions to the general spin-boson 
problem have applications in a huge range of systems, from 
quantum computing to physical chemistry and photosynthe¬ 
sis [41]. Various techniques and methods exist to numerically 
acquire such solutions, one of the most powerful of which is 
the hierarchy equations of motion [42,43]. Here we use those 
equations to model a two-level system coupled to a bosonic 
environment or reservoir. The general Hamiltonian is written 
as 


_ _ , V 

Hsb = —cr 2 + Act x + ^2 WkC^ak + EE ®+ a k J , 

k k 

(13) 

where A is the two-level system tunneling amplitude, and E 
is the two-level system splitting. The environment modes are 
described with creation (ajj and annihilation operators (a k ) 
with energy tu k , which couple to the system, described by the 
Pauli operators a z and o x , with strength Z k . By assuming 
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that the environment modes are well-described by a Drude- 
Lorentz spectral density, J(w) = 2 auj c , where a is the 

system-reservoir coupling strength and oj c is the bath cut-off 
frequency, we can exactly solve the dynamics of the two-level 
system (details can be found in Ref. [41-43]). We can then 
compare the non-Markovianity as measured via the temporal- 
S W to that given by the non-monotonic behavior of the entan¬ 
glement, as given by the concurrence [1], between the two- 
level system and an isolated ancilla [32]. One important dif¬ 
ference in the two approaches is that in the temporal-SW case 
there is no ancilla. In the ancilla case, the initial condition be¬ 
tween the system and ancilla is that of a maximally-entangled 
state; to mimic that in the temporal-SW case we assume the 
two-level system is initially in a maximally-mixed state. We 
then evolve the entire system-reservoir equations of motion, 
using parameters relevant to energy transfer in photosynthe¬ 
sis [41], and plot both measures in Fig. 2. 

For both measures we see similar behavior, particularly as a 
function of reservoir cut-off frequency and reservoir temper¬ 
ature. However, as a function of system-reservoir coupling, 
the entanglement measure has a larger window of detection. 
This may be attributed to the hierarchical relationship between 
EPR steering and entanglement. For example. Ref. [2] has 
shown that EPR steerable states are a superset of Bell non¬ 
local states, and a subset of entangled states. This hierar¬ 
chy links together these three different notions of quantum 
correlations. Therefore, the fact that the concurrence-based 
measure of non-Markovianity is more sensitive to the non- 
Markovianity than the temporal-SW measure seems linked, 
intuitively, to the notion that steering, in its EPR form, is a 
subset of entangled states. Also note that, the sharp features 
in both measures are typical, and arise because of the sudden 
vanishing and reappearance of both quantities in the temporal 
domain. Note that here, for consistency with Ref. [32], we 
plot TVtsw and Me using 


Mi 


r l f 

dfi[p(t)] 

>to 

dt 


dt + fi[p(tf)} - fi[p(t 0 )], 


(14) 


where for the temporal-SW measure i = TSW, the function 
fi[p(t)} is the temporal-SW at time f, while for the concur¬ 
rence measure, i = C, the function /, [p(t )] is the concurrence 
between system and ancilla at time t. This definition for the 
integral differs from Eq. (7) by a trivial factor of 1/2. 

Conclusions .— To summarize, we have discussed the con¬ 
cepts of “temporal” steering and how this can be quantified in 
a similar way to that of the original spatial EPR-steering. We 
further proved that the temporal steerable weight decreases 
monotonically under a CPT map and can be used as a measure 
of non-Markovianity, suggesting that both forms of steering 
can act as a quantum resource, similar to entanglement. Fi¬ 
nally we note that, in parallel, the temporal steerable weight 
has been recently implemented experimentally [44]. 
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FIG. 2. (Color online) The non-Markovianity measures, A/tsw 
(based on the temporal steerable weight) and Me (based on the entan¬ 
glement with an ancilla), as a function of system-reservoir coupling 
a, reservoir cut-off frequency cj c and temperature T, for a two-level 
system coupled to a bosonic reservoir with Drude-Lorentz spectral 
density. The system parameters are chosen to maximize the affect of 
the reservoir memory, with E = 0 and A = 100cm -1 . The mag¬ 
nitude of these parameters are typical for excitation energy transfer 
in photosynthesis [41], where the memory effect and structure of the 
environment is thought to play an important role. 
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In this supplementary material we give a few illustrative examples of the calculation of the temporal steerable 
weight and its application as a measure of strong non-Markovianity for some prototype models. 


HOW TO CALCULATE THE STEERABLE WEIGHT: 
A PEDAGOGICAL EXAMPLE 


Here we show explicitly how to calculate the steerable 
weight of Skrzypczyk et al. QJ in a simple example. Specifi¬ 
cally, we assume three types of measurements corresponding 
to the projections on the eigenstates of the Pauli operators: 


x = l+X+H-X-l, 

Y = \R)(R\-\L)(L\, 

z = |o)(0| 11)<11, (1) 

where |0> = \H), |1) = \V), |±> = (|0> ± |1»/V2, 
\R) = (|0) + z11))/x/2, and | L) = (|0) - i\l))/y/2, which 
can be interpreted as: horizontal, vertical, diagonal, antidi¬ 
agonal, right-circular, and left-circular polarization states for 
the optical polarization qubits, respectively. We can label the 
eigenstates of the Pauli operators together with their eigenval¬ 
ues as follows: |a;i) = |+) with x\ = +1, \x 2 ) = |—) with 
X 2 = —1, |yi) = |+) with yi = +1, ..., and | z 2 ) = |1) with 
^2 = -1. 

Then, possible unnormalized states of Bob <j a \ x (x = 
X , Y, Z) for a given two-qubit state p read 

^1 }l = °+i\x = Tr A [(|+}(+|(g)/)p], 

a a\l = a -Mx = Tr^fd—)(—|®7)p], 

°a\l = a +l\Y = 

a a\l = a -i \y = Tr A [(\L)(L\®I)p], 

a^ x =a +llz = Tr A mmm, 

a a\l = a -i\z = Tta[(|1)(1|®/)p], (2) 

where I is the single-qubit identity operator. A classical ran¬ 
dom variable held by Alice, 


A n — Uj ; — [(^z \X\xi ), (jji |Y|yj), (Zi \ Z\z \)\, (3) 


can take the following values: 

Ai = [-1,-1,-1], 
A3 = [—1, +1, —1], 
A 5 = [+ 1 ) — 1 ) —l]i 
A 7 = [+1, +1) —1], 


A2= [—1, —1, +1], 

A4= [“I, +lj +1], 

A6= [+1, —1, +1], 

A8= [+1, +1) +!]■ (4) 


The extremal deterministic single-party conditional probabil¬ 
ity distributions for Alice read 

[D Xl (+l\X), ...,D Xs (+1|X)] = [0,0,0,0,1,1,1,1], 

[D Xl (-1| A), ...,D Xs (-1|*)] = [1,1,1,1,0,0,0,0], 

[D Xl (-l\Z),...,D Xs (-l\Z)\ = [1,0,1,0,1,0,1,0]. (5) 

Let us denote an unsteerable assemblage as 

8 

^ = X D 8 ( a \ X ) a A = X Dx " ( a \ X ) a An • (6) 

A n= 1 

Then, we have 

( 1 ) US __ US 11 

a a\ x = a +i\x = + *>* + <*\T + <*>*> 

( 2 ) US _ US III 

a a\x = CT- 1 |X =Cr Ax +CTA 2 +Cr A3 +CTA 4 , 

( 3 ) US _ US 11 

a a\x = a +l\Y = ^3 +cta 4 +cr Xl + a Xs , 

( 4 ) US US _ 1 1 1 

a a\x = a -l\Y ~ a \i + 0A 2 + Ct a 5 + OA 6 , 

( 5 ) US _ US 11 

<|' = (T +1 | Z = (7a 2 +CTa 4 +CTa 6 +CTa 8 , 

( 6 ) US _ US 111 r 7 \ 

' = V-i\Z = a M + °A 3 + • (7) 

The steerable weight SW can be given as the solution of the 
following semidefinite program: Find 

8 

SW = 1 — maxTr^ cta„^ ( 8 ) 

n =1 

such that 

( a a\x ~ CT Si US ) > 0 and > 0 (9) 

for i = 1,2,..., 6 and n = 1,..., 8. By using a numerical 
package for convex optimization f2H4], one can implement 
this semidefinite program in a straightforward way. This is 
easily generalized to the temporal case by replacing the two- 
qubit measurements in Eq. ({2} with measurements on a single 
qubit, followed by evolution under the channel A. 

EXAMPLE 1: COHERENT RABI OSCILLATIONS OF 
A MARKOVIAN SYSTEM 


As a first simple example of the behavior of the temporal- 
SW under a Markovian dynamics, we consider a qubit that 
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FIG. 1. (Color online) The temporal steerable weight (temporal-SW) 
as a function of evolution time when a system is in (a) a Markovian 
environment (example 1) and (b) non-Markovian environment (ex¬ 
ample 2). (a) The temporal-SW when the system undergoes coher¬ 
ent Rabi oscillations and purely Markovian decay (example 1). The 
black dashed, red solid, and blue dotted curves represent the results 
of the decay rate 71/57 = 0, 1/6, and 1, respectively. The time 
t is in units of 1/gi, and Ti is set to 1. (b) The temporal-SW when 
the system interacts with a non-Markovian environment (example 2). 
The black dashed, red solid, and blue dotted curves represent the re¬ 
sults of the decay rate 72 /J = 0. 0.03, and 0.1, respectively. Here, 
the time t is in units of 1 / J. 

undergoes coherent Rabi oscillations and purely Markovian 
decay. The Hamiltonian of the system is 

H = hg 1 (a + + a-), (10) 

where hg\ is the coherent coupling strength between two 
eigenstates, |+) and |—), of the qubit, and = |+){—| 
and cr_ = |—)(+| can be considered the raising and lower¬ 
ing operators, respectively. A Markovian channel induces a 
dissipation rate 71 from |+) to |—). We assume that the ini¬ 
tial state, po in Fig. 1 of the main text, is a maximally-mixed 
state and then perform projective measurements M„\ x in three 
(or two) mutually-unbiased bases: X, Y, and Z (or X and 
Z). In Fig. 1(a), we plot the temporal-SW as a function of 
the evolution time t. We can see that the temporal-SW al¬ 
ways remains the maximal value of unity if there is no decay, 
while the temporal-SW decreases monotonically when 71 is 
non-zero, as expected; the dynamics of this system is Marko¬ 
vian. 


EXAMPLE 2: A SIMPLE NON-MARKOVIAN MODEL: 

A QUBIT COHERENTLY COUPLED TO ANOTHER QUBIT 

Our second example is that of a qubit coherently-coupled 
to another qubit. If we treat one qubit as the system and the 
other one as the environment (by tracing it out), we have a 



coupling strength g 

FIG. 2. (Color online) The degree of the non-Markovianity for a mul¬ 
timode reservoir with Lorentzian spectral density (example 3). The 
non-Markovianity A/Tsw, defined by the temporal steerable weight, 
as a function of the coupling strength g. Here, g is in units of spectral 
width oj w . 


very simple example of a non-Markovian environment. The 
total Hamiltonian of the system in the interaction picture is 

#int = hJ(a]_a 2 _ + cri<7+), (11) 

where a l + and u!_ are the raising and lowering operators of the 
zth qubit, and UJ is the coherent coupling between the system 
and the environment. We assume the system qubit is also sub¬ 
ject to an intrinsic decay with decay rate 72 . In Fig. 1(b), we 
plot the temporal-SW for various decay rates 72, after tracing 
out the effective environment-qubit. The initial condition of 
the system-qubit is that of a maximally-mixed state, while the 
environment-qubit is in its excited state. As seen in Fig. 1(b), 
there is a vanishing and a reappearance of the temporal-SW 
of the system qubit. Since we know that the temporal-SW 
should decrease monotonically under a Markovian dynam¬ 
ics, the oscillation of temporal-SW naturally shows that the 
qubit is undergoing non-Markovian evolution. This memory 
effect in this simple example is easy to understand in that in¬ 
formation regarding the state of the system-qubit flows to the 
environment-qubit and returns at a later time; one cannot as¬ 
sume that the evolution of the environment is not influenced 
by its history. 


EXAMPLE 3: A QUBIT COUPLED TO 
A NON-MARKOVIAN MULTIMODE RESERVOIR 

In general, the dissipation 7 rate in a Master equation de¬ 
scription of an open-quantum system can be time-dependent, 
i.e. 7 = 7 (f). If 7 (t) < 0, it indicates that information can 
flow back to the system and the system dynamics can be non- 
Markovian. To show that the temporal-SW is sensitive to this, 
we use the same example as in Breuer et al. Q], where a qubit 
is coupled to a reservoir with a Lorentzian spectral density. In 






3 


this case, the decay rate can be written as 


l(t) 


2 d 

G(t)di 


\G(t)\, 


( 12 ) 


where 


G{t) = e~^ t/2 




(13) 


with b = y/ui^ — 2 gu> w . Here, g denotes the coupling 
strength and uj w is the spectral width. We choose a mixed state 
as the initial state and plot the non-Markovianity A/tsw as a 
function of g/u> w in Fig. 2. Our results agree well with those 
in Ref. the non-Markovianity is zero when g/u> w < 0.5, 
and increases monotonically as a function of g/ui w . 
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